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Abstract 

This article provides novel analytical results for the Rice function, the incomplete Toronto function and the 
incomplete Lipschitz-Hankel Integrals. Firstly, upper and lower bounds are derived for the Rice function, Ie{k,x). 
Secondly, explicit expressions are derived for the incomplete Toronto function, TB{m,n,r), and the incomplete 
Lipschitz-Hankel Integrals of the modified Bessel function of the first kind, 7e,t,„(a,z), for the case that n is an 
odd multiple of 0.5 and m > n. By exploiting these expressions, tight upper and lower bounds are subsequently 
proposed for both TB{m,n,r) function and /e^,„(a,z) integrals. Importantly, all new representations are expressed 
in closed-form whilst the proposed bounds are shown to be rather tight. Based on these features, it is evident that 
the offered results can be utilized effectively in analytical studies related to wireless communications. Indicative 
applications include, among others, the performance evaluation of digital communications over fading channels and 
the information-theoretic analysis of multiple-input multiple-output systems. 

Index Terms 

Closed-form representations. Rice Je-function, Incomplete Toronto function. Incomplete Lipschitz-Hankel Integrals, 
Marcum Q-function, lower and upper bounds 

I. Introduction 

It is widely accepted that special functions constitute invaluable mathematical tools in the majority of fields in 
natural sciences and engineering. In the area of telecommunications, their utilization in various studies studies 
often allows the derivation of analytic expressions for important performance measures such as for example 
error probability and channel capacity. Furthermore, it has been shown that the computational realization special 
functions is not generally laborious since the most of them are included as built-in functions in popular scientific 
mathematical packages such as Maple, Matlab and Mathematica. To this effect, both the algebraic representation 
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and computation of any derived analytical expressions have been significantly simplified. 

Among others, the Rice /e-function, the incomplete Toronto function and the incomplete Lipschitz-Hankel 
integrals (ILHIs) appear in analytical solutions of numerous problems in telecommunications. They were all proposed 
several decades ago and they are denoted as Ie{k,x), TB{m,n,r) and Ze^,„(a, 2), respectively. 

The Rice /e-function was firstly proposed by S. O. Rice in [1] and has been largely exploited in the study of 
zero crossings, in the analysis of angle modulation systems, in radar pulse detection and in error rate analysis 
of differentially encoded systems, [2]— [5]. Its definition is typically given in integral form which involves an 
exponential term and a modified Bessel function of the first kind and zeroth order. Alternative representations 
include two infinite series and one closed-form expression - reported in [4] and [5], respectively. The former, are 
expressed in terms of the modified Struve function and the modified Bessel function of the first kind whereas the 
latter is expressed in terms of the Marcum Q-function, Qm{a,b), [6]— [10]. 

In the same context, the incomplete Toronto function constitutes a special case of the Toronto function, which 
was initially proposed by Hatley in [11]. It also includes as a special case the Marcum Q-function and has been 
used in studies related to statistics, signal detection and estimation, radar systems and error probability analysis, 
[12]— [14]. Its definition is also given in integral form while alternative representations include two infinite series, 
which were proposed in [15]. 

Finally, the ILHIs are a class of incomplete cyUndrical functions that have been largely encountered in analytical 
solutions of numerous problems in electromagnetics, [16]— [17] and the references therein. In communication theory, 
they have been sufficiently utiUzed in investigations associated with the error rate analysis of MIMO systems under 
imperfect channel state information (CSI) that employ adaptive modulation, transmit beamforming and maximal 
ratio combining (MRC), [18]. 

However, in spite of the evident importance of the Ie{k,x) and TB{m,n,r) functions and the Ze^,„(a, 2;) 
integrals, it is noted that they are all neither tabulated, nor included as built-in functions in the aforementioned popular 
software packages. As a consequence, they are rather inconvenient to handle both analytically and computationally. 

Motivated by these issues, this work is devoted in deriving novel representations for these functions and integrals. 
Specifically, upper and lower bounds to the Ie{k, x) function and explicit expressions and upper and lower bounds 
to the TB{m,n,r) function and /e^,„(a, 2;) integrals, are derived for the case that n -|- 0.5 e N and m > n. 
Notably, the offered results are expressed in closed-form and have a tractable algebraic form. As a result, they can 
be meaningfully utilized in various analytical studies associated to wireless communications such as the performance 
evaluation of digital connmunications over fading channels and the information-theoretic analysis of MIMO systems. 
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II. Definitions and Existing Representations 



A. The Rice le-function 



The Rice /e-function is defined as [4]— [5], 

Ie{k,x)= [ e-Ho{kt)dt, 
Jo 



0<k<l 



(1) 



where /o(.) is the modified Bessel function of the first kind and zero order, [6]— [7]. An equivalent integral 
representation to (1) was given in [4], namely, 



Ie{k, x) 



1 1 ^-ni-ro.o) 



\/l - fc2 TT Jo 1 - kcosO 
along with the following alternative series representations. 



de 



(2) 



Ie{k,x) 
1 



2\/l - A;2^ ;^n! 2"\/r^P^ 



(3) 



and 



leik, X) = xe-^^ £ f + 



xe 



n=0 



oc 



2fc ; r(n+|) 

'x{l-k^)Y In+i{kx) 
2k 



(4) 



„=ov - J r(n + i) 

The notations and r(.) denote the modified Struve function and the gamma function, respectively [6]— [7]. 



RecalUng [4], equation (3) converges relatively quickly for the case that a;\/l — k'^ is large and kx is small. On 
the contrary, equation (4) converges quickly when x\/l — is small and kx is large. Nevertheless, this way of 
computation of Ie{k,x) function is rather inefficient due to the following three facts: a) two relationships are 
required; b) the above series are relatively unstable due to their infinite form; c) the function is not built-in 
in widely used mathematical software packages. 

An adequate way of resolving this issue was reported in [5]. There, the Rice 7e-function is related to the Marcum 
Q-function of the first order, Qi (a, b) by the following relationships, 

1 



Ie{k, x) 



[2Q{a, 6) - e-^/o(A;x) - 1] 



and 



Ie{k,x) 



1 



[Q{a,b)-Q{b,a)] 



(5) 



(6) 



where, a = y/xyl + VT^W and b = ^/x\J \- \/\ — k'^. 
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B. The Incomplete Toronto Function 

The incomplete Toronto function is defined as, 

rs(m,n,r) ^2r"-'"+ie-''' / i"-"e-*'/„(2rt)dt (7) 







Importantly, for the special case that n = ^^^^h^, it can be equivalently expressed in terms of the Marcum Q-function 
by the following relationship, 



Tb (^m, ^^^y^, = 1 - Q™+i (rV2, bV2^ 



(8) 



Alternative representations to the Tsia, b, r) function, in the form of infinite series, were reported in [15], namely, 

n2a 2(n-a+l) , °° r^l'Y, 

TB{m,n,r) = j e (9) 

and 

rr f \ 2(n-a+i) -r'^r^''lia + k,B) 

TB{m,n,r)^r^ ^^-^^—^ (10) 

where the notations (a)/c and 7(0, x) denote the Pochhammer symbol and the lower incomplete gamma function, 
respectively [6] — [7]. Also, 



and 



m + 1 

Although equation (9) is exact, its algebraic representation is rather inconvenient both analytically and numerically. 
Equation (10) is significantly more tractable than (9), yet its infinite form eventually raises convergence and 
truncation issues. 

C. The Incomplete Lipschitz-Hankel Integrals 
The general ILHI is defined as, 

Ze„,,n{x;a)^ f x'"e-'^"Z„(a;)da; (11) 







where m,n,a,z may be complex [16]-[17]. The notation Zn{x) denotes one of the cylindrical functionj^ Jn{x), 
In{x), Yn(x), Kn{x), H}^{x) or H^{x), [6]-[7]. An alternative representation for the ILHIs of the first-kind modified 
Bessel functions, was recently reported in [18]. This representation is given in terms of the Marcum Q-function 
and is expressed as follows, 

m n+1 
i=0 j=0 

'Only the In{x) function is considered in the present analysis. 
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(12) 



where the set of coefficients Al^^ „(a) and Blj^j^{a) can be obtained recursively, [18]. As aforementioned, the above 
relationship was found useful in analytical investigations related to error rate of MIMO systems under imperfect 
chaimel state information (CSI). 



III. Performance Bounds to the Rice 7e-FUNCTiON 

This section is devoted to the derivation of upper and lower bounds for the Rice le-function. To this end, it is a 
critical to express Ie{k,x) function alternatively. 



A. An Alternative Representation to the Ie{k, x) Function 



Lemma 1. For x > and < k < 1, the following relationship holds, 



Ie{k,x) = l-e-''Io{kx)+k / e-Ui{kt)dt 

Jo 



(13) 



Proof. By integrating once equation (1) by part, it follows that. 



Ie{k, x) 



e~*dt 



loikt) 



e~*dt 



loikt) 



dt 



According to the basic principles of integration, the first integral yields straightforwardly, 

e-*dt = -e"^ 



/ 

Jo 



(14) 



(15) 



Subsequently, based on [6]— [7], the derivative of the modified Bessel function of the first kind and order n is 
re-written as, 

d k 

(16) 



d k 

— /„(A:a;) = - [/„_i(fca;) + In+i{kx)] 



which for n = reduces to 



d 
di 



Io{kt) = kh{kt) 



(17) 



Consequently, by setting x = t, substituting (15) and (17) into (14) and noticing that /„(0) = 1, equation (13) is 
deduced. ■ 
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B. An Upper Bound to the Ie{k, x) Function 



Theorem 1. For x > and < k < 1, the following inequality holds, 



Ie(k, x) <l-e ""loikx) + \l 



VT+k 



(18) 



Proof. It is recalled here that the modified Bessel function of the first kind is strictly decreasing with respect to its 
order, n. Therefore, for a > it immediately follows that In-a{x) > In{x). Thus, for n ± a ± ^ € N, one obtains 
the following inequaUty, 



Ie{k, x) <l- e ^/( 



{kx) + k 

Jo 



e~*Ii{kt)dt 



(19) 



Notably, for the special case that n is an odd multiple of 0.5, i.e. n + 0.5 e N, a closed-form representation for 
the In{x) function is given in [6, eg'. (8.467)], namely. 



fe=o ^/^fc!(r^-fc)!(2a;)'=+5 
Therefore, for n = 0, eq.(20) reduces to, 

h{kt) 



(20) 



(21) 



\/27rfci 

Evidently, the proof of the theorem is subject to evaluation of the integral in (19). To this end, by substituting (21) 
into (19), one obtains. 



^ e-Hi{kt)dt = 
which has the following closed-form solution, 



\/2nkt 



dt 



£e-*I^{kt)dt=^ 



rfiy^VT^) erfi^VTTk) 



VTTk 



(22) 



(23) 



where. 



A 2 /"^ 
erf{x) = -/= 

Jo 



e"* dt 



is the error function, [6]— [7]. Finally, by substituting (23) into (19), the proof is completed. 



Remark. The authors in [19] derived closed-form bounds to the Marcum Q-function. Thus, by making the necessary 
change of variables and make the according substitution in equations (5) and/or (6), an alternative expression to (18) 
can be deduced. However, this expression is significantly less compact and convenient than (18) both analytically 
and numerically. 
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C. A Lower Bound to the Ie{k, x) Function 



Theorem 2. For x > and < k < 1, the following inequality holds, 



leik, X) > 2Q(^ + ^)+2Q(6-a)-e-7o(M-l (^4) 

Vl — k"^ 



where Q{.) is the one dimensional Gaussian Q-function [6], 



Q{x) 4 -1= / e-Vdf 



Proof. According to the aforementioned monotonicity property of the /„ {x) function it follows that /| (a;) < 7i (a;). 
Therefore, by making the necessary change of variables and substituting in (13), one obtains. 



Jo 



Ie{k,x) > 1 - e-^'Ioikx) + k / e-H^{kt)dt (25) 



Importantly, a similar inequality may be also deduced by exploiting equations (5) and (6). To this end, it is firstly 
recalled that the Marcum Q-function is strictly increasing with respect to its order m. Based on this, it follows that 

Qi(a,6) > Qi(a,&) (26) 

Subsequently, by substituting (26) into (5), the following inequality is deduced, 

(27) 



Ie{k,x) > 2Qi{a,b)-e-'=Io{kx)-l 



Of note, the authors in [19] show that, 

Qi{a,b) = Q{b + a) + Q{b-a) (28) 
Therefore, by substituting (28) into (27), eq.(19) is obtained and thus the proof is completed. ■ 

Remark. A lower bound to Ie{k,x) could be theoretically derived by following the same methodology as in 
Theorem 1. Nevertheless, this approach ultimately renders a representation that is both complex and divergent. 



D. Numerical Results 



The behaviour and tightness of the derived bounds is illustrated in figures 1-4. In more details, figure 1 depicts 
the bounds in (18) and (24) for k = 0.5 along with results obtained from numerical integrations for comparative 
purposes. Evidently, the upper bound is tighter than the lower bound for small values of x. However, its tightness 
degrades as x increases while the lower bound becomes much tighter. According to figure 2, this appears to be also 
the case for a relatively small x and different values of k. 

Interestingly enough, for large values of x- typically x > 40- the lower bound in (24) is so tight that it eventually 
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becomes a highly accurate approximation to Ie{k,x). This is clearly illustrated in figure 3 where for a; = 80 and 
different k, the plotted curve provides an excellent match to the corresponding theoretical results. This is also 
justified by the overall small absolute relative error, Car =| Ie{k,x) — eg. (19) | /Ie{k,x), which is ear < 10~^ 
over almost the whole range of values of k. 



IV. A Closed-form Representation and Bounds to the Incomplete Toronto Function 

As mentioned in section /, the TB{m,n,r) is neither expressed in terms of other elementary and/or special 
functions, nor is a built-in function in popular mathematical software packages. Motivated by this, we derive a 
closed-form representation for the case that n is an odd multiple of 0.5. Subsequentiy, we exploit this representation 
to propose novel closed-form upper and lower bounds. 

A. A Closed-form Solution to the TB{m,n,r) Function 

Theorem 3. For m,r,B >0, m>n+^ and n -|- 0.5 e N, the following closed form relationship holds, 



TB{m,n,r)=Y.Y. 



(-1) 



m—k—l 



-'ot'o V^k\{n-k-^)\l\{L-iy. 



(29) 



/ + 1 



where L = m—n—k—h and r(a, x), 7(0, x) denote the upper and lower incomplete gamma functions, respectively. 



Proof. By setting in (20) x = 2rt, the /„ {x) function is re-written as. 



n— :7 

1 ^T- (n 



U2rt) = - E 



+ k+^)\r-''-i 
k\ (n-fc-i)!22fe+i 



(-l)'=e^'-* + (-l) 



n+ig-2rt 



n-l- - e N 



(30) 



By substituting into (7) and utilizing the basic identity (a ± 6)^ = ± 2ab + 6^, it follows that. 



TB{m,n,r) = ^ -r 



l\]„n—m-k- 



fe=0 



Jo Jo 



(31) 

/o Jo 
where L — m~n — k — j. Evidently, a closed-form solution to the above expression is subject to evaluation 
of the two involved integrals. To this end, with the aid of [21, eg. (1.3. 3. 18)], the above expression is re-written 
equivalently as follows, 

""^ ^ (n + k-\)\L\2-^<'r-^<'-' 



(-1) 



m—k- 



/ t'e-*'di- (-1)"*-'=-' / t^e-^'dt -{-!)'' t^e-*'dt+{-l)^ 

Jr J B+r J B—r J —i 



I B+r 

Finally, with the aid of [6, eg. (3. 381. 3)], eq.(29) is deduced and thus, the proof is completed. 



f'e"* dt 



(32) 
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B. Upper and lower bounds to the Tb {m, n, r) function 

With the aid of Theorem 3, explicit bounds to the incomplete Toronto function may be straightforwardly deduced. 

Corollary 1. For m,r,B >0 and n e N, m>n, the following inequality holds, 

TB{m,n,r) >TBi^,n + ^,r^ (33) 
where Tb [m,n + ^,r) is given in closed-form in eq.(29). 

Proof. The incomplete Toronto function is strictly decreasing with respect to n. Therefore, for an arbitrary real 
positive value a, it follows that TB{m, n + a,r) < TB{m, n, r). Hence, for the case that n e N and a = 0.5, the 
lower bound in (33) is straightforwardly deduced in closed-form. ■ 

Corollary 2. For m,r,B >0 and n gN, m>n, the following inequality holds, 

TB{m,n,r) <Tb (m,n-^,r^ (34) 
where Tb {m,n — ^,r) is given in closed-form in eq.(29). 

Proof. The proof follows immediately from Theorem 3 and Corollary 1. ■ 

C. Numerical Results 

The validity of the closed-form expression and the behaviour and tightness of the derived bounds are shown in 
figures 4 and 5 with respect to r. In the former, eq.(29) is depicted for n = 0.5 and m = 1.0 along with results 
obtained from numerical integrations for m = 1.0 and n = 0.4, n = 0.5 and n = 0.6. In the latter, eq.(29) is 
depicted for m = 3.0 and n = 2.5 along with numerical results for n = 2.4, n = 2.5 and 2.6. It is evident that (29) 
is in exact agreement with the numerical results while the overall tightness of the derived bounds is quite adequate 
over the whole range of values of r. 

v. a closed-form representation and performance bounds to the incomplete 

Lipschitz-Hankel integrals 

Likewise the TB{m, n, r) function, the ILHIs are neither tabulated, nor built-in in popular mathematical software 
packages. However, their algebraic form renders possible the derivation of a closed-form expression for the case 
that n is an odd multiple of 0.5. By exploiting this result, upper and lower bounds are deduced based on the same 
criteria as in the previous section. 
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A. A Closed-Form Solution to the Iem,n{z;a) Integrals 



Theorem 4. For m>n, m,r,B>0 and n + 0.5 e N, the following closed-form relationship holds, 



where P = m 



, 7(^,(a-l)^) , , 
' {a-l)P 



■1) 



„+i 7(J^.(a + l)^) " 
{a + ir . 



(35) 



Proof By substituting the function with its closed-form representation in (20), equation (11) can be equiv- 

alently re-written as, 



7e„,4^;«)-^^ x/^fc!(n-fc-i)! 



(-l)*^ / a;^e-"^e^da; + (-1) 



(36) 



Notably, the integrals in (36) clearly belong to the family of gamma special functions. Therefore, after some basic 
algebraic manipulation and with the aid of [6, eg'. (3.381.3)], eq.(35) is deduced. ■ 



Remark. The present analysis was hmited in the consideration of only the In{x) function in (11). However, similar 
results may be also derived analogously for the Bessel functions Jn{x), Yn{x), Kn{x) and the Hankel functions 

B. Upper and Lower Bounds to the Iem,n{z', a) Integrals 

Corollary 3. For m> n, m,r,B >0 and n the following inequality holds 

Im,n{z'^ «) > ^m,n+i {^'^ «) (37) 

where the I^.n+il^^^) is given in (35) in closed-form. 



Proof The Im,n{z]a) integrals are strictly decreasing with respect to n. Therefore, for an arbitrary real positive 
value a, it follows that Im+a,n{z] a) < Im,niz; a). Thus, for the case that n e N and a = 0.5, one obtains the the 
closed-form lower bound in (35). ■ 

Corollary 4. For m>n, m,r, B > and n e N, the following inequality holds 

Im,n{^'^ O) < Im,n-^{z; O) (38) 

where Im,n-^{z\o) is given in (33) in closed-form. 



Proof. The proof follows immediately from Theorem 4 and Corollary 3. ■ 
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C. Numerical Results 

The validity and behaviour of the offered results are exphcitly illustrated in figure 6. Specifically, one can observe 
the exactness of (35) and the rather tight performance of the performance bounds in all areas of values of their 
arguments. 

VI. Conclusion 

In this work, exphcit representations and performance bounds for the Rice /e-function, the incomplete Toronto 
function and the incomplete Lipschitz-Hankel integrals of the modified Bessel function of the first kind were 
derived. The offered results are novel and are all expressed in closed-form. This property is sufficiently advantageous 
since it renders them suitable for efficient utilization in various analytical apphcations in the wide area of digital 
connmunications. 
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Fig. 1. Behaviour of tlie bounds in equations (18) and (24) for k = 0.5 
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Fig. 3. Behaviour of the lower bound in (24) for x = AO 
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Fig. 4. Behaviour of tlie exact solution and the peiformance bounds to Tsim, n, r) for m = 1.0 and different n 
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Fig. 5. Behaviour of tlie exact solution and the peiformance bounds to Tsim, n, r) for m 



= 3.0 and different n 
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Fig. 6. Behaviour of the exact solution and the peiformance bounds to /em,n(z; a,) for different values of n and m 



